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Abstract: We compute the back-reaction of pairs of codimension-two branes within an ex- 
plicit flux-stabilized compactification, to trace how its properties depend on the parameters 
that define the brane-bulk couplings. Both brane tension and magnetic couplings to the sta- 
bilizing flux play an important role in the resulting dynamics, with the magnetic coupling 
allowing some of the flux to be localized on the branes (thus changing the flux-quantization 
conditions). We find that back-reaction lifts the classical flat directions of the bulk super- 
gravity, and we calculate both the scalar potential and changes to the extra-dimensional and 
on-brane geometries that result, as functions of the assumed brane couplings. When linearized 
about simple rugby-ball geometries the resulting solutions allow a systematic exploration of 
the system's response. Several of the systems we explore have remarkable properties. Among 
these are a propensity for the extra dimensions to stabilize at exponentially large sizes, pro- 
viding a mechanism for generating extremely large volumes. In some circumstances the 
brane-dilaton coupling allows the bulk dilaton to adjust to suppress the on-brane curvature 
parametrically below the change in brane tension, potentially providing a mechanism for re- 
ducing the vacuum energy. We explore the stability of this suppression to quantum effects 
in the case where their strength is controlled by the value of the field along the classical flat 
direction, and find it can (but need not) be stable. 
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1. Introduction 

Most of what is known about the physics of branes situated within extra dimensions either 
neglects their back-reaction onto their environment, or approximates the surrounding geom- 
etry as noncompact by ignoring the physics responsible for its stabilization at finite volume. 1 
Although these are often good approximations, there are also very interesting situations where 
they are not. 

A particularly interesting case where these effects cannot be neglected is when it is the 
back-reaction itself that stabilizes some of the extra-dimensional moduli. This case turns 

1 Randall-Sundrum models |ffl] are important exception to this statement, where back-reaction is incor- 
porated through the Israel junction conditions 0], but these are restricted to the limiting special case of 
codimension-one branes. 
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out to be important for compactifications whose extra-dimensional volume is very large, such 
as those arising within large- volume string vacua ||. In particular, the larger the extra 
dimensions the lower the string scale Mi, and once the string scale gets as low as the TeV 
scale — such as in supersymmetric extensions Q of ADD-type models || — supersymmetry 
becomes dominantly broken on the branes rather than by the fluxes in the bulk |7|, ]|| . In this 
case it is known that brane-induced corrections can dominate the leading classical predictions 
for the potential governing the lightest moduli || . 

The need to include back-reaction when computing the shape of the low-energy scalar 
potential is both a potential asset and a liability. The downside is the additional complexity 
required to properly incorporate both the extra-dimensional and brane dynamics within a 
controlled approximation. The upside is the potential for progress finding new mechanisms 
for understanding long-standing problems. Progress in particular on naturalness problems 
to do with the existence of light scalar masses and small vacuum energies, that hinge on 
understanding all contributions relevant to the low-energy scalar potential. And there are a 
variety of reasons for thinking that brane dynamics could be useful for understanding these 
problems Jl(], ||] . 

Six dimensional super gravities provide a fruitful place to explore these issues because they 
are complicated enough to exhibit many of the features of ten- or eleven-dimensional string 
vacua, yet they are simple enough often to allow explicit solutions and so more systematic 
exploration of the various configurations of physical interest. 6D gauged chiral supergravity 
[pH has proven particularly useful, providing early insights into chiral fermions and flux 



compactifications [|12|, 13 1. This has motivated finding a great many exact solutions to the 
classical field equations for this system, including a broad class of flux compactifications 
for which the two extra dimensions are a warped, squashed sphere with singularities at the 
positions of two positive-tension source 3-branes. These include solutions for which the on- 



brane geometry is flat [||, 14, [L5| (also known as 'rugby-ball' solutions), de Sitter/anti de 
Sitter like [16|, time dependent [jlTj or involves other bulk fields |l8| or additional branes fi~9[| . 

In this paper we explore brane back-reaction in this system by computing how the flat 
rugby-ball solutions respond to a general perturbation of the brane-bulk couplings. In par- 
ticular, we assume the perturbed brane-bulk couplings to be given by the leading terms in a 
derivative expansion of the brane action, 



s b = / lnu) + $ b *F), (l.i) 

where uj is the volume form for the space- filling 3-brane, and *T denotes the 6D Hodge 
dual of the background Maxwell flux, T MN (whose presence stabilizes some of what would 
otherwise be light bulk moduli, in the same way that 3-form fluxes stabilize some moduli 
in ten-dimensional flux compactifications J20||). The coefficient Tj, denotes the tension of the 
brane in question, which can be an arbitrary function of the bulk scalar dilaton, eft, appearing 



in 6D gauged chiral supergravity. ${, has a similar interpretation H 21 1 as an on-brane flux 
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and can sometimes compete with r& to play an important role in the low-energy energetics of 
the back-reaction. 

The bulk geometry that interpolates between a generic pair of source branes is known to 



be time dependent [17|, in much the same way that a random collection of mutually interacting 
electric charges is also not static. This is reflected by the generic absence of time- independent 
solutions once a brane-bulk system is perturbed. Unlike earlier stability analyses for these 
geometries [|22| , we do not deal with this by seeking the time-dependence of the solutions to 
the brane-perturbed bulk equations of motion. Rather, we instead couple an external current 
that stabilizes this time-dependence in order to study the energetics of the potential energy 
that drives it. In practice, at low energies this current need only couple to the massless 
Kaluza-Klein (KK) 'breathing' mode of the leading-order extra-dimensional geometry, since 
this is a flat direction in field space along which the time dependence dominantly lies. 

In this way we find the response of the on- and off-brane geometries as a function of the 
perturbing brane couplings, as well as the shape of the scalar potential that stabilizes and gives 
a mass to the low-energy breathing mode, for general choices for the brane coupling functions 
Tb and <3?fc. We find instances where the breathing mode is stabilized by the interaction of the 
branes on the bulk, as well where it instead runs away to infinity (which, perhaps 

surprisingly, includes the simplest case where both and <3?6 are independent of the 6D bulk 
dilaton, </>). 

When restricted to the special cases for which our results duplicate earlier calculations, 
we fully reproduce earlier expressions. But our systematic survey of perturbed solutions also 
reveals some new ones with surprising properties. These include (see §4 for a more detailed 
summary): 

• Solutions whose extra-dimensional volumes stabilize at values that exponentiate any 
moderately large hierarchies among the brane-bulk couplings, naturally giving enor- 
mously large volumes; 

• Solutions whose on-brane geometry can be parametrically small compared with the 
largest energy scales that appear in the brane-bulk couplings (though, alas, not yet 
small enough to describe the observed Dark Energy density); 

• Solutions for which the value of the breathing mode along the low-energy flat direction 
defines the strength of both brane and bulk loop corrections, and for which this ensures 
that the above two properties can be stable against quantum effects; 

• Models for which the brane-bulk couplings can have the form required to profit from a 
'chameleon' mechanism 



Our presentation is organized as follows: The next section, §2, describes the linearized 
solutions to the bulk field equations, and how the integration constants in these solutions are 
determined by matching to the functions Tj, and <3?6 that define the codimension-2 bulk-brane 
interactions. These are then used to provide explicit expressions for the extra-dimensional 
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and on-brane geometries as functions of these brane properties. The results of the full 6D 
calculation are compared with the effective 4D picture that captures the low-energy limit, 
since the scalar potential in this effective theory provides an efficient way to understand the 
implications of brane dynamics on low-energy properties. This section closely follows the 
logic of ref. |^], which performs a similar calculation in the non-supersymmetric case. 

§3 then uses the general results of §2 to explore the implications of several simple illustra- 
tive choices for the coupling functions r& and A particularly simple toy model — for which 
*l2b T b K X/b^*£> c*- f° r small r\ — is also examined, that exhibits modulus stabilization at 
exponentially large volume and parametric suppression of the low-energy on-brane curvature 
(or vacuum energy) . Finally, this section estimates the effects of brane and bulk loops for the 
toy model, and argues that the exponentially large volume, and the small on-brane vacuum 
energy (and scalar masses) can be technically natural. 

Our conclusions are summarized in S4. 



2. The bulk-brane system 

This section defines the system of interest. The fields of interest are part of the bosonic sector 
of chiral gauged supergravity in six dimensions [jllj, for which we follow the implications of 
coupling to nonsuper symmetric branes. In particular we follow the metric, g MN ', a bulk 
Maxwell gauge potential, A M , whose presence helps stabilize the bulk geometry; and the 6D 
scalar dilaton, cj). 

2.1 Field equations and background solutions 

We first describe the bulk equations of motion and brane boundary conditions, followed by a 
simple class of rugby-ball solutions near which general solutions are sought. 

Bulk equations 

The bosonic action in the bulk is 2 

S huik = -J dW^S ^9 MN {n MN + d M ^d N ct>)+-^T MN T MN + ^eA , (2.1) 

where the two dimensionful constants are the gauge coupling, g R , for a specific U R (1) sym- 
metry of the supersymmetry algebra, and the 6D gravitational constant, k. One of these sets 
the overall scale of the bulk physics, leaving the dimensionless combination g R / k as a free 
parameter. Here T = dA denotes the gauge potential's field strength. 

The equations of motion from this action are the (trace reversed) Einstein equations 

Tl MN + d M 4>d N (f> + K 2 e-+F UP F N P - e-^ PQ T PQ - e^j g MN = , (2.2) 

2 We use a 'mostly plus' metric and Weinberg's curvature conventions j24|] (that differ from those of MTW 
] only by an overall sign in the definition of the Riemann tensor). 
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the Maxwell equation 

V M {e~*F MN ) = 0, (2.3) 



and the dilation equation 



n4>-^ie't>+ l ^e-'t>T MN T MN = 0. (2.4) 
k z 4 



Since these field equations are invariant under the transformations 

9m n -»• C,9mn and e - * -»• Ce"^ , (2.5) 

with A — > A, any nonsingular solution is always part of a one-parameter family of solutions 
that are exactly degenerate (within the classical approximation). 

Symmetry ansatz 

In what follows we restrict attention to solutions that have maximal symmetry in the four on- 
brane directions, and axial symmetry in the two extra dimensions. This assumption restricts 
us to solutions involving at most two source branes. The corresponding ansatze for the metric 
and Maxwell field are 

ds 2 = dp 2 + e 2B d6 2 + e 2W g^dx^dx" and A = A d9, (2.6) 

where g^uix) is a maximally symmetric metric, and all of the functions W, B, <p and Aq 
depending only on p. The corresponding Maxwell field strength is T^q = A' , where primes 
denote differentiation with respect to the coordinate p. 

Subject to this ansatz the bulk field equations reduce to 



(e-B + 4W e -<p A , d y = {Aq 



( e B+4W 0')' _ f M e <t> _ I K 2 Q 2 e 4> e ~%w\ e B+4W = q 



(0 



W" + (W') 2 ] +B" + {B') 2 + {^') 2 + - A K 2 Q 2 e't > e~* w + 9 -^e't > = (pp) (2.7) 

B " + ( B ')2 + m > B > + 3 2 Q 2 ^-SW 5| ^ = 

4 K 

i e-^i? + W + 4(H^') 2 + WB' - \ k 2 Q 2 e^e- 8W + 4 e* = . 
4 4 k 

The first of these can be integrated once exactly, introducing an integration constant, Q, 
labeling the bulk flux, 

F pe = A' e = Qe*e B - AW . (2.8) 
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Rugby ball solutions 

In the special case that the dilaton is constant, <j) = (fo, these equations have a simple solution 
with extra dimensions having the shape of a rugby ball, sourced by two branes ||: 



ds 2 = e"^ 



dp 2 + a 2 L 2 sin 2 (^)dO 2 



+ -q^dx^dx" 



Toe = ^^ o/2 = Qe^ o/2 aL 



sm 



(2.9) 



where g uu = r} av denotes the usual flat metric of Minkowski space. The extra-dimensional met- 
ric becomes singular at the brane positions, p N = and p s = ttL, which are (^-independent 
because of the coordinate rescaling 

The geometry generically has a conical singularity at these points, characterized by the 
defect angle 5 = 2tt{\ — a). In the special case a = 1 the extra-dimensional geometry is a 
sphere, corresponding to the supersymmetric Salam-Sezgin solution [O]. The deficit angle 



can be related to the common tension, T, of the two source branes by [26] 



a 



K 2 T 

2vr 



(2.10) 



The equations of motion impose two relations amongst the integration constants, requir- 



ing 



-2 



K 2 Q 2 



K 

and k 2 Q 2 L 2 = 1 



(dilaton equation) 
(Einstein equation) 



(2.11) 



Additionally, flux quantization due to the spherical topology of the extra dimensions implies 



- = 2aL 2 Q = — , 
9 9r 



(2.12) 



where g is the gauge coupling of the background Maxwell field and n is an integer. The 
couplings g and g R are in general different because the background Maxwell field need not 
be the one that gauges the U R (1) symmetry. This last condition determines the deficit angle, 
a, and thereby constrains the tension of the source branes. As is elaborated in more detail 
below, a minor modification of these solutions allows the source branes themselves to carry 



some of the total flux, ^branes; in which case eq. (|2.12|) generalizes to 



n 
9 



Q! ^branes 

9r 2vr 



1 

9r 



k 2 T 
~2n 



+ 



branes 

2vr 



1 

9r 



2tt 



branes 



(2.13) 



where the last equality uses eqs. ( pTH ), which imply Q = 2g R / 'k 2 . This can be regarded 
as allowing the tension in these solutions to be arbitrary, provided the on-brane flux is also 
dialed, $branes (T), to satisfy eq. Q2.13f) . 
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For fixed brane flux the above construction describes only a one-parameter family of 
solutions, labeled by ipQ. This one-parameter degeneracy is the one required by the scale 
invariance, eq. (2.5), of the classical field equations. Because of the overall factor of e^^ 
in the extra-dimensional metric, eq. ( |2.9D , the proper distance between the two branes is 
Ap = e -v> °/ 2 7rL and the volume of the extra dimensions is 



V 2 = 4iraL 2 e 



(2.14) 



Our interest in what follows is in how this flat direction gets lifted by dilaton couplings 
to the branes. Its connection to the extra-dimensional volume makes this also a stabilization 
mechanism for the size of the extra dimensions; a codimension-2 generalization of the better- 



known Goldberger-Wise stabilization mechanism for codimension-1 branes [27| within RS 
models. 



Brane matching conditions 

We take the brane-bulk coupling to be defined by the following lowest-derivative action, 
including both a (^-dependent tension and a ^-dependent coupling to the Maxwell field ||: 

1 



Si 



branes 



b=N,S 



b=N,S ' 



-9 e 



AW 



(2.15) 



where the coupling functions r b and <& b can depend on all of 4>, W and gee without breaking 
the condition of maximal symmetry in the on-brane directions. Because of the explicit factor 
of extracted from the Maxwell coupling, these interactions also do not break the bulk 
scaling symmetry, eq. (|2.5| ), only when both T b and & b are (^-independent. Our conventions 
are such that e pe = \/y/g~2 = e~~ B transforms as a tensor, rather than a tensor density, in 
the two transverse dimensions. The parameter r b has the physical interpretation of being 
the tension of the brane, and (as is shown below) the parameter & b similarly denotes the 
magnetic charge (or flux) carried by the source branes. 

The presence of such brane couplings imposes a set of boundary conditions on the deriva- 
tives of the bulk fields in the near-brane limits, 3 given by 4 

* 2 T 6 



B 1 1 



with Th ■- 



and 



1 



Pb 



in, 



Pb 



U b with 
- T b + 3U b 



2vr 

2 <m 



(2.16) 



3 These matching conditions can be derived [g9| from codimension-1 microscopic models for codimension- 
2 branes. 

4 Notice that we normalize the quantities % and lAb without including the factor of e 4W used in this reference. 
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where, as before, primes denote differentiation with respect to p. is defined in terms of if, 
and &b as the total lagrangian density of the source, 



T b :=n- Qbe-fe-"?,*. (2.17) 

As shown in Appendix [A|, the corresponding boundary condition for the Maxwell field 
implies that the integral of T p e to obtain Ag(p) in a coordinate patch containing each source 
brane gives 

Ag(p) = — + Q ^ dp e ^ B ~ m Northern hemisphere 
27r Jpn 
rp 

+ Q dpe^ +B ~ 4W Southern hemisphere , (2.18) 

J PS 



2vr 

where <3?b := lim^^ $>b[4>(p)] — appropriately renormalized |28j — and the signs are dictated 
by the observation that increasing p points away from (towards) the North (South) pole, 
together with the requirement that the two patches share the same orientation. Requiring 
these to differ by a gauge transformation, g~ 1 doil, on regions of overlap implies the flux- 
quantization condition 

n_£tot +Q f Ps dpe 4>+B-AW j (219) 



9 2vr 



PN 



which identifies $tot = as the part of the total magnetic flux carried by the branes ||. 

2.2 Perturbations 

In this section we use the previous discussion to analyze how couplings to the brane lift the 
flat direction associated with the scaling symmetry of the bulk theory, and so to see how the 
scalar zero mode, ipo, becomes stabilized at a specific value, (fo = tp*. Our discussion closely 



follows the discussion of the nonsuper symmetric system in ref. [21]. 

It is instructive to contrast how this stabilization differs from the nonsupersymmetric 
system. To this end recall how the stabilization occurs in detail, from the point of view of six 
dimensions. Given two branes, we seek the bulk configuration satisfying the field equations 
that interpolates between the boundary conditions that each brane specifies. Specializing to 
solutions that are both axially symmetric in the transverse directions and maximally sym- 
metric in the on-brane dimensions requires seeking bulk profiles that depend only on p. 

What is important is that the brane boundary conditions only specify the derivatives 
of the fields near the branes, and not the values of the fields themselves there. Once the 
derivatives of the fields are specified at one brane, the values of the fields at the same brane 
can be adjusted to try to ensure that the derivatives take the values required by the other 
brane at the other brane's position. It is in this way that the stabilized value, cpo = <p+, is 
obtained if the brane actions break the classical bulk scaling symmetry. 

This argument shows that a classical solution satisfying all of the boundary conditions is 
in general impossible given an arbitrary choice for tpQ. From the low-energy 4D perspective 
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the absence of a solution when ipQ ^ tp* corresponds to the absence of a static solution for 
a value of (po that is not an extremal of the low-energy effective potential, V^ s ((po) ^ 0. It 
can still be possible to map out the shape of the scalar potential for generic (po, however, 
provided we turn on an external current, J, coupled to <po that is designed to ensure that 
(po is a stationary point of the potential, including the current. The shape of the effective 
potential can be computed by seeing precisely how much current is required as a function 
of ipo. In what follows we define the current coupling by adding the following term to the 



Sj = - d«xy/=g J , (2.20) 



action 5 



where J is a constant (since our goal is only to couple a current to the would-be zero mode, 

In this kind of construction the stabilized value, 99*, corresponds to the choice for which no 
external current is necessary, </(<£>*) = 0. An important difference between the supersymmetric 
system of interest here and the nonsupersymmetric one studied in ref. is that in the 
supersymmetric case it can (but need not) happen that there is no value of (po for which 
J((po) = 0. As we shall see, from the 4D point of view this corresponds to an effective 
potential that is a pure runaway, for which V^ s ((po) only vanishes as (po — > ±00. 

Linearized equations 

Our goal is to solve the above field equations by linearizing them about a rugby-ball solution. 
This amounts to assuming that the (^-dependent contribution to 77, is small relative to the 
tension that is responsible for the rugby-ball geometry itself: 

T b = t + Snitf,) and $ 6 = $ + S$ b {<t>) , (2.21) 

with the background deficit angle sourced by T = r — Q The linearized equations of motion 
including the current term — derived in Appendix (|C]) — are given below. All background 
(rugby-ball) quantities are denoted by a subscript 0, and perturbations are universally denoted 
by 5: so Q = Qo + SQ etc. Since we ignore all second-order quantities we may write 5Q/ Q ~ 
SQ/Qo and so can use either of these quantities interchangeably. Also, since Wo = for the 
rugby balls, W = 5W . 

To linear order the Maxwell field strength becomes 

T pe = QaLe^' 2 sin (J£\ (l + ^ + SB - ASW^j , (2.22) 



the on-brane curvature is 
R = -Ae w 



25W +±cot( I) d s 5W + dlSW 



L 2 L \L p p 



+ (2 . 23) 



L 2 \ Q 



As is shown in Appendix j3|, most of the low-energy physics of interest is insensitive to the detailed form 
of the current to which we couple, so long as it has a good overlap with the would-be zero mode. 
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and the remaining linearized field equations become 



sin 

. o 

sin 



dp(SB) 



1 

Z 2 



sin 



(t) 



1 

L2 



Q 1 S1 " 



3 / SQ 



l(^L)_QSW + K 2 JL 2 e-w 



4 

— — cot 
1j 



dp5W 



1 



and dfSW = - cot [ ^- ) d h 5W . 



P 



(2.24) 



Finally, the linearized flux quantization condition can be expressed as 



6Q 
Q 



l 



2.3 Linearized solutions 



dp sin 



41 (4-5VF - <5£ - <50) - ^ f 5$^ + 5<S> S ) . (2.25) 
L / 47TO V / 



The strategy is to construct the general solution to these linearized equations, and then to 
use the brane matching conditions to eliminate the resulting integration constants in terms 
of brane properties. To simplify expressions it is convenient to define the dimensionless 
coordinate x := p/L = (p/^e^^ ^ 2 , keeping in mind that its implicit dependence on <fo 
brings this dependence to any bulk fields that depend on x. We have some freedom in how 
to group the perturbations; which we employ (without loss of generality) to simplify the 
linearized flux-quantization condition as much as possible. 

First, we solve the equation for the warp factor, W, which has the general solution 



5W(x) = W + Wi cos x , 



(2.26) 



where Wo and W\ are integration constants, of which Wo = may be ensured by rescaling 
the on-brane coordinates, x^. 

With this solution, the equation to be solved for the dilaton becomes 



sin x d x (5(p) 



5Q S 

4W\ cosx — | siu .)■ 



which integrates to give 

5<j) (x) = 5tpo + Vi hi 



cosx 



smx 



— 2W\ cos x + 



6Q 
Q 



In I s'mxl 



(2.27) 



(2.28) 



Here Stpo an d fi are integration constants, of which 5<po = can be ensured without loss of 
generality by absorbing it into the otherwise arbitrary background value, (po ■ 
Finally, the equation of motion for 5B becomes 



d x [sin 2 x d x {8B)} 



sin 2 x 



-<pt In 



1 — COS X 


SQ ( 


sinz 


Q \ 



+ In I sinxl 



+12VFicosx - K 2 JL 2 e 



(2.29) 



-fo 
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This integrates to 



SB = 5B + B\ cot x - 4Wi cos x - (piM.2(x) 



(2.30) 



+■ 



5Q 
Q 



- x cot x — H2(x) 



+ - (k 2 JLV w ) xcot; 



where 5Bq and B\ are integration constants. Of these, B\ is pure gauge in that it can be 
changed arbitrarily by reparameterizing the coordinate p. We fix this coordinate freedom by 
defining p = to be the position of the 'north' brane, which requires e B — > as p — > 0; 
ensuring B\ = 0. The functions M.2 and H2 appearing here are defined by 



Mi(x) :-- 
M 2 (x) :-- 



dy sin yln 
M x {y) 



1 — cos y 



smy 



dy 



■ 2 ' 
sin y 



(2.31) 



and 



%\(x) := / dy sin 2 y In | siny| 
Jo 



n 2 (x) 



dy 



sin 2 y 



(2.32) 



For later convenience when discussing flux quantization it is useful to absorb parts of 
these integrals into the definition of 5Bq, by writing 



8B = 5B - 4Wi cos x + <pi 



1 



^-M 2 (x) 



+ 



5Q 
Q 



- x cot x — H2(x) 



+- (k 2 JLV^ ) (x cot x + 1) , 



with the number M defined by 



M 



dx sinxA^aO 



(2.33) 



(2.34) 



Numerically this evaluates to the value 6 JVl = — 1, which we use throughout what follows. 
Flux quantization 



Using the above expressions in the linearized flux quantization condition, eq. ( 2.25 ), gives 

K 2 Q 



5Q 
Q 



dx sinx 45W - 5B 



5® N + <5$c 



-SBo + 



3 [5Q\ k 2 Q 



Q 



5® N + 5®c 



(2.35) 



"Maple 11, 10 digit precision, see Appendix (|d]) 
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which uses the integral 



n 



dx sinx% 2 (x) 



-0.613706 ~ln4- 2. 



(2.36) 



and the last approximate equality is a numerical inference. 7 The absence of ip\ on the right- 
hand side of eq. ( p.35[) is a consequence of the definition of 5Bq used in eq. (fE33l) . Solving 
this for 5Bq gives 

_1 f5Q\ k 2 Q i 
4 V Q J 4vra 



<5<3?,v + 5<I>< 



(2.37) 



Finally, the linearized field equations return the following on-brane curvature, 
~25W 1 ( p 



R 



_4 e vo 

/.-' 

2e^° /<5Q 



I 2 \S 



L 



2k 2 J. 



L 



+ ^ cot \ ^- \ d a SW + d\ 5W 



L 2 \ Q 



(2.38) 



Notice that all of the ^-dependence cancels in this expression (as must happen given our 
assumption of maximal symmetry), leaving a result that is determined purely by the change 
of bulk Maxwell flux and the applied current. 

2.4 Physical interpretation and renormalization 

The above solutions are described by four physical integration constants, which we can take 
to be ipo, (pi, Wi and 5Q/Q. These can be traded for four physical properties of the bulk 
and on-brane geometries. 

Wi can be taken to be the difference between the value of the warping (which controls 
the gravitational redshift) between the two branes, which is given by 

5W N - 5W S = 2 Wi . (2.39) 

Similarly, to linear order the near-brane geometry as x = p/L — > is governed by 

1 + 8B(x) 

S /AO\ _ .1 

(2.40) 



e B ~ e _¥,o/2 aX sinx 
~ ap 



1 + S B -4Wi-^ + ^(^)+ K*JL 2 r-* + Oijr < 



which corresponds to a conical singularity (since the e B vanishes linearly with p), having 
defect angle a N = a + 5 a N with 



5 a> 



a 



SB - Wi - f + + . 2 JL 2 e~^ 



(2.41) 



7 Mathematica 7, with thanks to Ben Jackel. 
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By contrast, as x = p/ L — > tt we have 



(2.42) 



which uses Hii^ — e) = (1 — Ln4)(7r/4e) + O(e ). In particular, this shows that e B does 
not vanish at p = ttL. Instead e B vanishes at p = tt(L + 5L), indicating a change in proper 
distance between the branes: p s — p N = ir(L + 5L)e~ V0 / 2 . The amount of the change is 
obtained by comparing eq. ( fOg) to the Taylor expansion of e B about its new zero, giving 



5L 
T 



(2.43) 



The singularity at the 'south' brane is also conical (at linear order), with defect angle 
^iven by a s = a + 5 a s with 



(2.44) 



One might imagine that a further observable could be the difference between the value of 
the dilaton field at the two branes, 4> N — 4> s , since this governs the relative strength of some 
of the bulk couplings to each brane (such as the strength of the bulk Maxwell couplings to 
brane- localized charged particles). However a subtlety arises in this case because the profile 
4>{p) diverges in the limit that p — >■ p N and p — >• p s . For this reason we defer a discussion of 
this quantity to the next section, which deals with renormalizing these divergences. 



Brane matching and renormalization 

Ultimately the bulk integration constants should be related to physical properties of the 
branes that are the source of the bulk geometry; this is where the brane matching conditions 
play a role. In order to perform this matching we must specify a functional form for the 
brane tensions, 77,, and fluxes, <&b- We take both of these to be smooth functions of 4>, and in 
many (but not all) examples we imagine these functions to be extremized at 4> = that is, 
(dn/dcp)^^ = 0. 

The problem in practice with matching is that the argument of Tb and $6 is (fib '■= <ft(Pb), 
but the profile 4>(p) given in eq. ( 2.28 ) diverges as p — > pb- For instance, for x = p/L = e -C 1 
and x = tt — e we have 



and 



• (x = tt — e) = (fo + 2W\ + ipi In 



e 






+ ( 


2 


2 






+ ( 


£ 





5Q 
Q 

5Q 
Q 



ln\e\ +0{e). 



(2.45) 



This divergence is dealt with by renormalizing the parameters that define the functional 
form of Tb and and in particular those parameters that determine the values Vb- It can be 
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absorbed in the definitions of i>b by defining renormalized quantities, vj,: 



v N 



v s 



^ In (e) - ^ In (e/2) 
^ln(e) + ^ln( £ /2) 



(2.46) 



where the first expression is relevant at x = (the north brane positon), and the second one 
at x = ir (the south brane). With these definitions, 



lim 

e-»0 



and lim 

e-H)L 



e) - Ujv 
7r - e) - v s 



ip - 2Wi - v N 
ipo + 2Wi - u s , 



(2.47) 



and so t\,{4> — Vb) = Tb(tpo i 2Wi — Ui>) and so on. This is a useful redefinition because our 
interest really is in the value at which the zero mode, <po, gets stabilized, rather than on the 
value of (f) itself at the brane position. And this is finite as e — > with renormalized quantities 
(like Vb) fixed. 

With this construction in mind, there are four independent matching conditions: 
\e B d p B] =l-^k and \e B d p B] T = -1 + ^ , 

L P J p=0 27T IP 1 P =ttL 2lT 



(2.48) 



where, as before, Tb = Tb — Q <&b e~ iW ( p *>^ . The difference in signs between north and south 
brane arises because increasing p points away from the north brane but towards the south 
brane. 

Specialized to the dilaton profile, eq. ( |2,28] ), the first two of the above conditions become 



a \<pi + 



while the latter two evaluate to 



6Q 
Q 
6Q 

Q 



a 



a 



Q 

ipi 3 / 5Q 



f d5T N 
2vr V dcj) 
k 2 /dST 



SB + K 2 JL 2 e~^ 



(2.49) 



lib, - — + -1-^] +SB ( , + k-JL-<- 



-ipo 



K 2 

2ir 



S ) 



(2.50) 



where the change in brane action from the background value, T, is 5Tb := Tb — T = 5rb — 
Q5® b + Q®[4:5W(pb) - 5Q/Q]. However, the terms involving SW and 5Q/Q in K 2 5T b may 
be dropped in the matching conditions because their contributions are suppressed by an 
additional factor of k 2 QQ/2tt relative to the leading contributions. Hence, from here on we 
take 5T b ~ Sr b - Q5<S> b - 
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Eliminating 5Bq using eq. ( 2.37 ), and solving the above matching conditions gives 



6Q 
Q 

Wi 



Aira 

K 2 

Aira - 
lQira 

K 2 



5T' N + 8T' S 



STL - 6T' a 



1 



5T N + - STL 



5T S + I ST' 



Ana 



6T N + -6T'-Q5$ f 



+ [6T S + -5T' S -Q6$ S 



, (2.51) 



where 5T^ denotes ddT^/dcf). These expressions allow the elimination of the three integration 
constants (921, W\ and 5Q/Q) and the current, J, to be completely expressed in terms of 
brane properties and tpQ. 

In particular, the condition J = is satisfied when = (p*, where 



ST N + 1-ST'- Q5$ N ) + ( 6T S + ^ST'-Q6$ 



0. 



(2.52) 



¥>()=¥>* 



This expression determines the stabilized value, = 92*, as a function of the properties of 
the branes. 

On-brane curvature 

Finally, the curvature in the on-brane directions, regarded as a function of ipo, becomes 



7raL 2 e^° 



1 



R(<po) = -{5T' N + ST' S 



2ira 



K 2 L 2 Je-^) 



5T N + 5T S + -{5T' N + 5T' S ) - Q{5$ N + 5® s ) 



(2.53) 



Of particular interest is this result specialized to the value, <po = ip+, that solves the field 
equations in the absence of the current J (if such a value exists - more about this below). 
The curvature evaluated at this value is the curvature predicted by the field equations for the 
brane geometry, and eq. (2.52) allows it to be written in two equivalent ways: 



R 



5T N + ST S - Q{5$ N + 5$ g ) 



<P0=<P* 



(2.54) 



The second of these agrees precisely with the corresponding expression obtained in the non- 
supersymmetric case studied in ref. |21]. However this is not also equal to the first equality of 
eq. ( p. 54 ), because in the nonsuper symmetric case eq. ( |2.52| ) no longer holds, being instead 
replaced by 5T' N + 8T' S = 0. 

From the point of view of a brane observer this must agree with the (maximally sym- 
metric) curvature that is predicted by the 4D Einstein equations given a 4D vacuum energy, 
0eS- 



R = -4k 4 g eS , 



(2.55) 
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where K4 is the 4D gravitational coupling, given in terms of the 6D coupling, k, by 

1 4vraL 2 e-^ 



(2.56) 



Comparison gives 



R ( iraL 2 



l(8T' N + 5T>) . (2.57) 



2 



Notice that this agrees (to linear order) with the more general exact classical result obtained 
in eq. (3.81) of ref. 

Q eS = -Y J {u b + l -T^ , (2.58) 

given that U b vanishes to linear order. 

2.5 The low-energy 4D effective theory 

This section constructs the effective 4D theory that reproduces the low-energy dynamics of 
(po and the 4D metric predicted by the full 6D theory. We do so at the purely classical level, 
working perturbatively about a rugby ball solution, as above. 

General form 

In this section the two fields of interest in the low-energy theory are the 4D metric, g^ u , 
describing the massless KK graviton, and a 4D scalar, 8 ip, describing the low-energy would- 
be zero mode, (po, associated with the scaling symmetry of the bulk field equations. (We 
ignore here any other low-energy fields, such as other 4D scalars or 4D gauge fields coming 
from A M or the metric.) 

The most general possible local 4D effective theory describing the interactions of ip and 
9fiu, up to the two-derivative level, is 



'eff 



J d 4 x^j [t u [f^R^ + H^d^d^] + V JF (<p)+jk(<p)} , (2.59) 



where /, h, Vjp and k are all functions to be determined, and j denotes a low-energy current 
that is included to explore the shape of these functions (in precisely the same manner as 
J was included in the 6D theory). Our task is to identify these functions by matching the 
predictions of this theory with the low-energy predictions of the full 6D system. 

The functions /, h and k differ from V w in that they already receive their leading contri- 
butions when the two source branes are described by their background tensions, T; without 



8 We use ip to denote the 4D field in the effective theory, to distinguish it from the (closely related) parameter 
ifo appearing in the 6D solutions. 
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the symmetry-breaking, (^-dependent contributions 5Tb(4>). These leading contributions can 
be obtained by simple dimensional reduction, which predicts 



and k(ip) = e~ v if we define joc4vraL 2 J. (2.60) 

To the same approximation the (Jordan frame) potential vanishes, Vjp(ip) = 0, since the 
background branes do not break the classical bulk scaling symmetry. 

Low-energy matching conditions 

The goal is to determine how these quantities are perturbed by the addition of ^-dependence 
to the brane action, (<p). Our main focus is on the contribution to Vjp, since (unlike for 
the other functions) for Vjp this is the dominant contribution. We use the prediction for the 
low-energy scalar curvature, R, as a function of ip — i.e. eq. ( |2.53| ) — as our means for doing 
so. 

To make the comparison we compute R in the low-energy effective theory, assuming 
a maximally symmetric geometry. Defining for notational convenience l/ftf := e^* / n\ = 
4naL 2 / 1 k 2 , the metric and scalar equations of motion are 

4k| 



'-§-je^ + yj>)=0. (2.61) 



Eliminating the current between these two equations gives the following expression for R as 
a function of ip, 

R ,„ d 



4% = V ' + r ' = ^6j(* V ')' (2 - 62) 



To obtain V w we regard eq. ( 2.62 ) as a differential equation to be integrated with respect 
to (p, using eq. ( 2.53Q to evaluate the left-hand side as an explicit function of ip. The integral 
yields 

V JF {p) = i e-f J dp (sT N + 5T S — Q 5$ N - Q 5<S> S + ^ST' N + | ST'A (2.63) 
f 5T N + 6Ts) + l - e -f J dp e* (-Q6$ N - Q 5$ s + ~6T' N + i 6T' S J . 



The integration constant, C, implicit in this integration contributes an amount C e - ^ to V JF , 
with C fixed by matching to the 6D theory at a specific value of (p. A convenient place for 
doing so is the vacuum configuration (if this exists), ip = defined by j(ip+) = 0, for which 
a prediction — eq. ( 2.52j ) — is known in the 6D theory. 



-17- 



Specifically, solving eqs. fl2,61|) for j(ip) gives 



je-' p = -[VU<p)+2V JF (<p)\, (2.64) 

and so (p+ satisfies 

V' JF (^) + 2V JF {^) = 0. (2.65) 

This has a simple interpretation in the Einstein frame, which is defined by rescaling = 
e,vf—'P*)g go that the 4D action has a canonical Einstein-Hilbert term 

S cS = ~ J d A x^g~ | gT [iV + 5 3^ d u <p] + Vxr (<p) + J ef } , (2.66) 

with j m := je~ 2tp * and 

V^-e^'^V^). (2.67) 
Clearly cp* therefore satisfies Vj£p(y>*) = 0, as might have been expected. Imposing V' JF +2V JF = 



when ip = satisfies eq. (|2.52| ) then gives 



■<5T 6 '(^) + Q^ 6 (^ 



8Tl(<po) - Q5$ b (<Po) 



(2.68) 



Given the Einstein-frame potential, classical vacuum energy is 
as found earlier (using eq. ( 2.52D ). The scalar mass similarly is 



(2.69) 



5 L 



(2.70) 



Similarly, chasing through the earlier expressions for the shape of the bulk geometry gives 



5a.b 
a 
6L 
T 



5W N - SW* 



2ira 
3k 2 



lQna 



K 

K 2 r 



8na 



1 



6T S + -ST')-[ 5T N + - 5T' h 



1 



87ra 



(2.71) 
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Notice in particular that no warping arises unless the two branes carry different amounts of 
localized flux. This is by contrast with the nonsuper symmetric case [21], for which net warping 
always accompanies a tension difference for the two source branes. But in the super symmetric 
case the flux quantization condition does not allow such a tension difference without some of 
the flux being forced onto the branes. 



3. Illustrative examples 



The previous formulae with which the previous section closed represent the main results of 
this paper. We now explore their consequences through a number of illustrative special choices 
for the (/^-dependence of the tensions on each brane. 



3.1 Dilaton-independent tensions and fluxes 

Consider first the simplest example: where both quantities 5t{, and 5<&b are independent of 
(p. In this case the condition, J(y*) = 0, defining c/?* degenerates to 

J2( ST b~Q^b) = ^2{Sn-2Q6^ = 0, (3.1) 

b b 

so two situations need to be distinguished. Either a solution to the condition J = exists - 
which requires 5Tb = 2Q J2b ^b — or it does not. Consider each of these in turn. 



When J = has solutions 

If the constant quantities 5Tb and 5&b satisfy the condition J2b^ T b = 2QX)&<^&> then max- 
imally symmetric solutions to the 6D field equations exist for any value of (p+. Because no 
particular value of <po is selected, this shows that the flat direction that ipo parameterizes is 
not lifted. This is consistent with the observation that the brane action scales the same way 
as does the bulk action — and so does not break the bulk scaling symmetry — in the special 
case where 5Tb and 5&b are both (^-independent. 

In this case formulae ( [2.69D , ( [2.70D and (|2.71|) degenerate to £ e ff = rrvL = 5L/ L = 0, while 
eqs. ( |2.71|) reveal 5a.b = K 2 5Tb/2ira, as usual, and 5W N — 5W S = k 2 (5T n — 5T s )/87ra. The 
new perturbed solution in this case is a special instance of the general solution to the full 
nonlinear equations [14, [l5|, 16], all of which are known for the symmetries of interest to us. 
In particular, the assumption of constant brane action, 5T^ = 0, is known to be sufficient to 
ensure g c g = 0, while 5T N ^ 5T S induces warping. As initially argued in ||, it is the freedom 
to have nonzero on-brane flux, <i>b, that prevents the flux quantization condition from being 
an obstruction to reaching these solutions as perturbations to the initial rugby ball (as one 



might naively have thought ] ft2|] , if eq. ( 2.12 ) were read as forbidding the possibility of having 
perturbations to a&, and hence also to T&). 
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When J / cannot be avoided 



The perturbative solution found here also allows an exploration of what happens in the more 
general situation where the fluxes and tensions are not related to one another by Ylb ^ T b = 
2QJ2b^^b- In this case there is no choice for (po = p* that can ensure J((f+) = 0, implying 
that no solution exists at all to the linearized field equations, subject to the assumed axial 
symmetry and on-brane maximal symmetry. In this case studies of linearized stability |2^| 
and exact time-dependent solutions [[H]] suggest that the relevant solutions are necessarily 
time-dependent . 

We now show how this expectation for time-dependence can be made more precise in 
the present context, since J / implies the absence of a stationary point to the (Einstein- 
frame) scalar potential, Vef ((f), for any finite value of (p. To show this we must reconsider 
the expression derived above for V JF , but without using the condition V^p(^) = to fix 
integration constants. For (^-independent 8rb and <53>& expression ( |2.63 ) for V JF becomes 



V JF (<p) = \ e~* f dp e vo (ST N + 5T S - Q 5<Z> N - Q6$ 3 + ^ST' N + | 6T^j (3.2) 
= lY.{ 6T b-Q 6 ®b)+Ce-*, 

b 

where C is the integration constant in question. 

A natural choice for C is to demand that V JF remain bounded as e v — > 0, since this 
corresponds to the weak-coupling limit for which both <f) and A M do not strongly self-interact 
in the bulk. More precisely, inspection of the 6D action, eq. ( |2 . 1| ) , shows that the bulk 
scalar potential vanishes in this limit, allowing the constant part of <fi to be absorbed into 
the definition A M ■= C^Am- This argues that Vjp should not become unbounded in this 
limit, leading to the requirement C = 0. 

With this choice the Einstein-frame scalar potential becomes 



V&te) oc e 2 ^ Y, ( 6T b ~ 2 5 ®b) = e 2lp fa - 2Q 5$ 6 ) , (3.3) 

b b 

which describes a runaway to ip — > ±oo — whose sign depends on the sign of J2b(^b — Q fi^b) ■ 
The absence of a solution here to = for any finite value of ip is what underlies the need 
for a time-dependent solution from the perspective of the low-energy 4D observer. 



3.2 Dilaton-brane couplings, vacuum energy and volume stabilization 

The next paragraphs explore some of the implications of nontrivial brane-dilaton couplings. 
Of particular interest is how the bulk and brane geometries depend on the choices made for 
these couplings. We start with the case where 5rb and 5&b vary only weakly with <p, and 
move on to more strongly varying examples. 
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Linear dilaton-dependence 

Consider therefore the simple situation where both brane tensions and fluxes are linear in ip, 
with 

n = + r bl ip and $ b = $ b0 + $ 61 (p , (3.4) 

with Tftj and <&bi constant. Since many — though not all — physical quantities depend only 
on the average brane action and flux, T e g := \^ b T b and § c r := ^Ylb^b, it is useful to 
phrase our assumptions in terms of these, which have the form 



where 



T eS (ip) =T + T 1 ip and $ cff (</?) = $ + $i¥> > 
Tl '-= \ Yl ( Tbi -Q^bi) and $i := - <S> bi 



(3.5) 



(3.6) 



b=N,S b=N,S 

We describe the resulting geometry as a perturbation about a rugby ball solution, char- 
acterized by a background tension, T, and brane flux, $(T), related by the background 
flux-quantization condition, eq. ( |2,13 ), 



r- Q$ 



2vr 

K 2 



1 



9 



(3.7) 



With this choice, the condition J = defining becomes 

= 5T cS (^) - Q^eff(^) + ^ 5Tls(<P*) 
= (T -Q$o)-(T-Q$) + (T 1 -Q$ 1 )^ + 



whose solution, 



1 



(r -Q$o)-(T-Q$) + 



(3.8) 
(3.9) 



G*i-Ti 

in this case exists so long as Ti 7^ Q^i- 

We remark in passing that the assumed linear coupling does not preclude the existence 
of a vacuum configuration, ip = 93*, contrary to what happens for the nonsuper symmetric 
situation described in ref. [21]. What is different in the nonsupersymmetric case is that ip* 
satisfies J2b ^bi^*) = — rather than ^2 b (5T b + | <5T b ' — Q 5&b)^ =ip =0 — which has no 
solutions if Xl&^Uv 9 ) i s a linear function of (p. 

The Jordan- frame scalar potential, eq. ( 2.6§| ), in the 4D effective theory then takes the 
simple form 



Vnriip) = Q$i + (Ti - G*i)(p - </>*) - (Ti + G* x )e-^-^ 
and so the Einstein-frame potential becomes 



(3.10) 



V EF {tp)= Q$i + (Ti-Q $i)(</>-</>*) 



= 2(>-</?*) 



(Ti + Q$i)e 



(3.11) 
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Requiring the potential to be bounded from below implies T\ > Q&\. Notice that at <p = ip* 
this satisfies V^t/?*) = automatically (by construction), and ¥^(9?*) = —T\ there — which 
agrees with — \ J2b ^Kv 9 *) = ~^es( i P*)j as ^ must. The physical parameters computed 
from Vef using eqs. fl2~69p and ( ^70|) in this case therefore are 



£cff 



-T\ and m't 



K 2 

f (3Tl 



Q$i) 



(3.12) 



while the extra-dimensional response of eqs. ( 2.71 ) becomes 



5a.], 
a 
SL 
L 



2ira 



and 



3K 2 g eff 
8ira 



W„-W» 



3k 2 T x 
8ira 



2ira 



n(tp*) - Q$b(<p* 



8ira 



(3.13) 

(3.14) 
the condition 



For potentials that are bounded from below — i.e. those with T\ > 
T\ > suffices to ensure m 2 > (and q c h < 0). 

Three important properties of these expressions bear special emphasis. 

First, g e ff quite generally depends on the background quantities T and only through the 
combination T — Q<3? whose value is constrained by flux quantization, eq. (3/7). Consequently 
g e Q does not change at all as T is varied, because flux quantization demands <3? must also 
be adjusted in a way that precisely compensates. Any value of T is equally good, and what 
counts for physical predictions is only the extent to which the values T e g (<£>*) — Q<& e ff(¥>*) 
differ from the flux-constrained background combination, T— Q <F This property also remains 
true for the more complicated examples discussed below. 

Second, it is relatively easy to arrange ip+ ~ —50 using only a mild hierarchy of parameters 



on the branes. But eq. ( 2.14 ) then ensures that the volume of the extra dimensions, V2 = 
^TraL 2 e~ 2v * , is exponentially large compared with the intrinsic scales on the branes and in 
the bulk. 

Third, what is most striking about this example is that the size of g e g and m 2 is com- 
pletely independent of T, Q<1>, To and <2<3?o- In this way this example captures part of 
the more general magic of codimension-2 constructions; they can admit classical solutions 
- like the rugby ball itself — for which large tensions coexist with flat (or weakly curved) 
on-brane geometries. Why is the result independent of the (^-independent part of T e g and 
<£ c fj? Quite generally, we know from eq. ( |2.69| ) that g e g = Y^bi^bif-k) ~ Q^bi^*)] = 
2 [<5T c ff (</?*) — Q^csi^P*)], and so (apart for the special case where <5T e g- cancels QQ c s) the 
reason g e g can be small even when To — Q<3?o is large is because the condition J = drives 
(p+ out to such large values that the terms To — Q <3?o and (T\ — Q 3>i)</?* mostly cancel in g e g. 

One is drawn from this last observation to try to identify how robust this property is, 
both to the shape assumed for Srb(p) and to the size of radiative corrections. 
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Power-law brane actions 

In the previous example | ip* | becomes very big if To and Q $o are much larger in magnitude 
than are T\ and Q $1, and so the assumption that Sti, is linear in ip typically cannot be justified 
simply as the first term in a Taylor expansion. It is useful therefore to examine slightly more 
complicated functional forms for 5Tb(<p) and 5&b{<p) in order to probe the robustness of the 
previous example. 

Let us consider branes of the general form 

T b = T b0 + T bri ip r > and $ 6 = $ 60 + *6j7 <P v 1 (3-15) 

again with constant tu and <£>{,,. The effective brane action and flux, defined as before by 
T efi ■= \ ^2b T b and $b ■■= \ Y,b then give 

T eS (cp)=T + T rj <p r > and $ efi (p) = $ + % <P V , (3.16) 

where 

Ti := ~ ^ -Q$bi) and ^ := - <S> bi . (3.17) 

b=N,S b=N,S 

As before we perturb about a rugby ball solution with background tension, T, and brane flux, 
related by the background flux-quantization condition, eq. (|3.7| ), and so 

6T eS = (T -T) + T v <p r > and Q5$ eS = Q($ - *) + G^V* ■ (3.18) 

We find ip+ by using the condition J(</?*) = 0, or 5T c fj — Q5<3? c fT + l^eff = ^' which m 
the present case gives 

(To-Q$ )-(r-Q$) + ^r 1 [(^-Q^)^ + |^] =0. (3.19) 

Approximate solutions are possible when \ (T V — Q^ v )(p^\ 3> \r]T ri /2\, in which case the field 
stabilizes approximately at 

v ' = {t^qwT- (320) 

where D is defined by 

D = Q($ - *) - (To - T) := -(5T - Q5<S> ) . (3.21) 

This is a real solution if the signs of D and — Q ^ are the same. For 77 > it is also large 
- and so justifies a posteriori making the large- y?* approximation — if \D\ S> \T„ — Q&n\- 
With this solution we find the low-energy cosmological constant is 



£eff = -tfZj, (p*) = -vT^ 1 = ~ ' IV " ^-f^ • (3-22) 



r]T v D ( T v -Q $ r ^ 1/ri 
T v -Q$ v \ D 

This reduces to the cases previously considered in the special cases r/ = and 77 = 1. Writing 
£> e ff = ~{ r ]T) /ip i ,)/(l — Q^ r] /T ri ) shows this result is generically suppressed relative to D within 
the approximations used, since these include |</?*| 3> 1. Because tp cx [D/(T V — Q& v )] 1 / ri , with 
all other things equal this suppression becomes stronger for smaller 77 > 0. 
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Exponential branes 

As our final example, consider several commonly occurring cases where the brane action 
depends exponentially on 92. A simple case of this type is when the entire tension and flux — 
i.e. both background and perturbation — involve a common exponential, Ti{ip) = t^q + A^ e aLp 
and QMv) = Q$bo + B b e a v. 

In this case the average brane action and flux, T e g := ~ Ylb^b and 3> e ff : = \ Ylb nave 
the form 

T eS ( l p)=T + Ae a * and Q$ c ff(¥>) = Q<& + B e aip , 



where 



To = - (m ~ Q <S>ob) , $o = 2 Yl $of " 

b b 

A=l^2^- B b) and B = \Y. B *>- 



(3.23) 



(3.24) 



For instance, given the explicit factor of in the definition of the brane-flux coupling, 
eq. ( 2.15 ), the special case a = 1 and Ab = <&bo = (and so <3?o = A + B = 0) corresponds to 
having brane actions that do not directly couple to the bulk scalar <f>. 

As before we perturb about a rugby ball solution with background tension, T, and brane 
flux, <3>, related by the background flux-quantization condition, eq. Q2.13 ), 



ng R 
K* \ g 

The perturbations about this background become 

ST cS = (T - T) + A e atp and Q 8<S> e $ = Q($ - &) + B e a ^ 



(3.25) 



(3.26) 



The condition J(</?*) = defining ip+ as usual is ST e g (</>*) + |<5Tg ff (y>*) — Q£<J> e ff = 0, 
which in this case becomes 



A 1 + 



(3.27) 



where -D := Q(&o — $) — (To — T). This has solutions if the sign of both sides is the same. 

(3.28) 



The low-energy Jordan- frame potential, eq. fl2.68|) , then is 



V JF (<p) = de«P + C 2 e~* + C 3 
1 



with C\ 



Co 



a + 1 



3o . , 

1 + y ) A - B 



(2 + o)Ci 



_(A + \aA-BY' a (A + \aA-B) 
and C3 = — D , 



l+l/a 



D l+l/a 
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leading to a similar expression for the Einstein-frame potential, Vgp = V JF e 2 ^ v *\ At (p = ip+ 
the cosmological constant becomes 



QeS 



(Tp — T) — Q($p - $) 

1 + (a/2) - 
(Tp - Q$ ) - (2tt/ K 2 )(1 - ng R /g) 
l + {a/2)-B/A 



(3.29) 



where the last equality uses the value of T — Q $ dictated by flux quantization. The scalar 
mass at the extremum is similarly 



a {6 + 3a-2B/A\ 



10 



2 + „-2 B M r'F- T °>- g (*- # °> 



a ^6 + 3a-2£/4^ 2 



io V 2 + 0-257.4 / 



27T , 

(T -C*o)--a(l 



(3.30) 



Eq. ( 3.29j ) identifies three potential mechanisms for suppressing 



QcS- 



1. The first is if a —> 0, in which case the brane actions become (/^-independent and 
92* recedes to infinity. This is the suppression already encountered in the examples 
presented above. 

2. The second is if the (^-independent parts, To and $0) are related to one another in the 
same way as flux quantization imposes on the background values, T and <J>. (For the 
special case where the bulk flux is chosen to lie in U R (1) direction (so g = g R ) and n = 1, 
the background lies in the same flux category as does the super symmetric Salam-Sezgin 
solution |12|], for which T - Q $ = 0.) 



3. Finally, the third potential suppression occurs even if T — Q<I> / 0, provided |<B/»4| 3> 1- 
That is, if 5<3? e g dominates 5T e Q then this only affects the value of leaving £ c g, as 
always, of order <5T e ' ff ((/?*). 

3.3 Quantum corrections and technical naturalness 

All of the calculations of brane-bulk interactions provided in previous sections are performed 
purely within the classical approximation. As such they leave open the question of how 
robust their conclusions are to modification by quantum corrections. And since the streets 
are littered with classical examples having small vacuum energies, a proper treatment of 
quantum corrections is the crucial to any credible mechanism for understanding the small 
size of the observed vacuum energy. 

In this section, we take a small step towards filling in this missing step, more in the spirit 
of indicating a promising line of inquiry than in providing a polished example. Our interest is 
in quantifying the stability of both the size of the low-energy cosmological constant, g e Q , and 
the size of the bulk volume, V2, (when this is large compared with more microscopic scales). 
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The starting point is an enunciation of the essence of the problem: once parameters 
are chosen to ensure a large value for Q e g and/or V2, are these choices stable against the 
renormalization that results when heavy fields are integrated out? In extra-dimensional brane 
models this question necessarily has two parts, to do with integrating out heavy field on the 
brane and in the bulk. 

We here use one of the previously discussed examples as a toy model for estimating the 
size of quantum corrections. We choose a model that has both has an exponentially large 
volume and a small 4D on-brane curvature — i.e. vacuum energy 9 — and estimate the 
size of quantum corrections. The main idea behind this model is that it is the bulk field <j) 
itself that counts both bulk and brane loops, with weak coupling corresponding to <fi being 
large and negative. The influence of loop effects is then simply incorporated by tracking the 
(^-dependence of the quantum-corrected (1PI) action, for which the above arguments about 
brane-bulk back-reaction can be applied. 

A toy model 

The theory of interest is one of the 'power-law' models described earlier. For our starting 
point we take a background rugby-ball geometry whose background tension and flux satisfy 
the flux-quantization condition, eq. (|3.7|), 



K 2 



1 



ng R 

9 



(3.31) 



with the right-hand-side being small enough to allow semiclassical reasoning, but not tuned 
to be inordinately small. Such geometries have flat on-brane directions, and as above we seek 
to see how brane-bulk interactions modify this, including loops. 

For the perturbations to this geometry we choose the classical brane-bulk Lagrangian to 
have the power-law form, 

5T cS = T* (-4>) r ' , (3.32) 

with < T) < 1 (and the smaller rj is, the larger the suppression in g e s)- Here T* is a function 
of all of the on-brane degrees of freedom, tp, such as 

T* = // + 8^ + M V + A V 4 + • • • , (3.33) 

which defines the scale [i. In the vacuum ij) = and so = [i 4 . 

A hierarchy is dialled in by choosing the scale fi in T* to be small compared with the 
typical brane scale, M: i.e. fj 2 <C M 2 <1/k. Notice that taking —<p 3> 1 does not affect the 
mass of the ip particle (or other brane particles in general) at the classical level, because 
appears only as an overall factor in the brane action. The goal is to show that the energy scales 
set by q 1 Jq and V 2 ^ 2 can be hierarchically different from M, and that this can be protected 



9 In the model £> e ff is small inasmuch as it is parametrically suppressed relative to other scales, though not 
small enough numerically to describe the observed Dark Energy. 
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from quantum effects. Since V2 turns out to depend exponentially on T*, a relatively small 
hierarchy between T* and M 4 suffices to generate very large volumes. 

The classical part of the story is worked out above, with (choosing <£ e ff = 0) eq. ( [3.20D 
implying 



T - Q$ 



i/t) 



2vr 



ng R 



i/v 



from which eq. ( 2.14| ) gives the bulk volume, 



V2 = 47raL 2 e v * ~ AiraL 2 exp 



2vr 



1/7?- 



(3.34) 



(3.35) 



Eq. ( |3.22D similarly gives the on-brane vacuum energy as 



QeS - V T *(-^ 



rj(T - Q$) ^ 2vrr/ 



ng R 



2vr 



ng R 



-1/7, 



(3.36) 



revealing a power-law suppression of g e g relative to 1/k 2 , whose strength improves the smaller 
77 gets. (e.g. for 7/ = 1 this gives £ e ff ~ while T) = \ implies £ e ff k2 /" 8 > and so on.) 

Loop corrections 

We now argue that the choice /i <C M underlying the classical hierarchy is technically natural. 
We do so using the observation that it is the expectation of the bulk zero-mode, 99, itself 
that controls the size of these loops, so loop corrections can be incorporated into the above 
argument by making a modified choice for the (^-dependence of T c q . 

Brane loops: 

To see why this is so, imagine first computing quantum corrections involving loops of the 
on-brane field, ip. When computing these loops it is useful first to adopt a canonical normal- 
ization, ip ific := (—if) 11 / 2 ip, after which the strength of the self-coupling becomes revealed 
to be A c V>c with 

Xr = ~ ~TTT • (3.37) 



(-<p)V 

More generally, because (—f) 11 pre-multiplies the entire brane action, for the purposes 
of power-counting brane perturbation theory it plays the role of 1/h. This ensures that 
each additional loop is parametrically suppressed by an additional factor of (—(p)~ v , with 
dimensions made up using the typical brane scale, M. In particular, integrating out a heavy 
field of mass M should give a Wilson action (or, alternatively a calculation of the 'quantum' 
1PI brane action) of the form 



eff 



)" + Ti + 



T 2 



+ 



(3.38) 



and so on. 
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Here the T n generically depend on the brane fields, 10 much as did T*. The point of 
quantum hierarchy problems is that — on dimensional 11 grounds — each of the T n is gener- 
ically of order M 4 , rather than the smaller fj, 4 . The question is whether this ruins the above 
conclusions about the size of and so also of V2 and g e ff ■ 

Now comes the main point. Because each loop correction is suppressed by an additional 
factor of (—(J))"* 1 , none of them has the same (^-dependence as does T*. In particular, none of 
them require the vacuum value of T* also to be of order M 4 instead of fi 4 . Better yet, having 
T\ ~ M 4 3> X* ~ fi 4 keeps 99* stabilized at large negative values, enforcing the dominance of 
the leading, classical, approximation. 

To see this in detail we repeat the above classical calculation of the potential for the bulk 
modulus if using the loop corrected action, eq. Q3.38 ), rather than the classical expression, 



eqs. ( 3.32| ) and ( 3.33 ). For simplicity we take <3? e fT = also at the loop level, though none 
of our conclusions would change if we were to assume a similar loop expansion for the brane 
flux, 

- *l + 7^ + • • • . (3-39) 

with the dimensions of the $ n again set by the large mass, M, circulating in the loops. 

To one-loop order, the new terms in r e g- are independent of 4>, and so their modifications 
to the brane-bulk back-reaction are encompassed by the analysis given in the previous section. 
Assuming — cp* S> 1 this gives 

f D\ lh r,D . AS 

-9?* ~ — and g eS ~ , (3.40) 

\ T *J {-<P*} 

where D = (T — Ti) — Q($ - <f> x ) ~ (T - Q $) - T x is of order the larger of M 4 or T - Q $ oc 
2tt/k 2 . Since both of these scales are much larger than T+ oc /u 4 , the classical assumption 
that —if* is large (and all that comes with it) is not undermined by one-loop corrections. We 
assume here that D and T* share the same sign. 

The size of the two-loop correction can be similarly estimated. The equation that deter- 
mines (p± is, at two- loop order 

- D+r *<-*K i+ &) + ^( i -&) =!0 - (3 - 4l) 

Using \D\ ~ IT2I ~ 0{M 4 ) and (—tp*) 71 ~ M 4 //i 4 3> 1 shows the 2-loop term to represent a 
small correction to the value predicted for (— cp*). 

The two-loop contribution to g e ff comes in two parts. The first of these is through the 
change in though because \S(p^/ip^\ <C 1 this contribution is subdominant to the value for 
£> e ff already computed at one loop. In particular, it doesn't ruin the suppression of g e & by the 
factor 1 /y*. 



10 Although these would be local for the Wilson action, they need not be for the 1PI action [ p3[ . 
11 These are often stated to be of order the 'cutoff' scale, but we make the more conservative statement that 
they scale with the physical mass M because cutoffs generically cancel in all physical quantities |3^]. 
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The second type of 2-loop contribution to g e s comes from the fact that T e ' ff now includes 
a new term of the form 



-77-1 



yT 2 
-<p+ 



1 



(3.42) 



which should be compared to the original contribution, ~ rjD/cp*. Clearly the new term is 
subdominant if \D\ ~ | X2 1 S> T*. And so it goes for higher loops, each of which is suppressed by 
an additional factor of (— tp*)~ n — fi 4 /M 4 , by virtue of the stability of the initial stabilization 
at large values of — (p*. 

Notice that the brane loops also don't have a large relative impact on the light scalar 
mass, which eq. ( |2.7C| ) gives to be 



2 ^4 

m v - -5- 



- ST^fa) + 35T> S (^) - QS&afa 



(3.43) 



ij-i 



\5M*J 

( 3r]T A( „ ^-i 



3r]D 

^ 37?r 2 r^ 



rjM 4 f 



M p 2 vm; 



(up to one loop) 



(two-loop term) 



with the final estimates using T 2 ~ D ~ M 4 and T* 
loops: 

The previous loop estimates are restricted purely to brane loops because they rely on the 
assumed form of the brane-bulk coupling. But the back-reaction of the brane loops onto the 
bulk is computed classically (for the bulk theory) just as before. How big might be quantum 
corrections in the bulk sector? 

An estimate for the size of bulk loops can be made in a manner very similar to the one just 
used for brane loops, because e 2 * is the loop-counting parameter for the bulk 6D supergravity. 
The simplest way to see this is to re-scale the 6D metric according to g MN — > g MN := e~^g MN , 
in terms of which the action of eq. fl2.1|) becomes 



'bulk 



d b x 



9 e 
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4--n Mp n NQ T r T 4- 

+ ^ y y •> mn-> pq 1 



2yi 



(3.44) 

where C, is a constant. This shows that for bulk perturbation theory it is the constant value 
of e 2 ^ that plays the role of of h. (The same also remains true once the action's fermion 
terms are included ||11|| .) Each loop involving bulk fields therefore contributes an amount 
proportional to an additional power of e 2< ^, which is small when (j> is large and negative (also 
the regime of weak brane coupling). 

Now imagine integrating out fields in the bulk that are heavy relative to the KK scale. 
Loops of these fields potentially modify both the brane and bulk actions by new local in- 



teractions [35, 36]. The loop-generated couplings arising in this way cannot depend on ip in 



-29- 



the same way as does the classical action, (2.1), again indicating that these classical terms 
are not themselves renormalized. Loop-generated terms necessarily involve new interactions 
whose (^-dependence can be organized into a series in powers of e 2ip [B7[ . 

This leads one to expect that each bulk loop is exponentially suppressed, by powers of 
e 2ip * oc 1/V| when — (p+ 3> 1. In particular, these corrections to physical properties would 
therefore be expected to be sub-dominant to the brane loops considered above. 



4. Conclusions 



This paper computes the back-reaction of a pair of 4D codimension-two branes onto the 6D 
geometry that they source, within a framework of flux compactification that allows a complete 
calculation of modulus stabilization. Although performed with a particular (gauged, chiral 
[0]) 6D supergravity, the mechanisms exposed by our calculations rely only on broad features 
(like the presence of a dilaton and scale invariance of the classical equations) shared by a wide 
variety of higher-dimensional supergravities. This leads us to expect them to have a wider 
domain of validity than the particular 6D system studied here. 

The main calculational assumptions are these: (i) we assume all energy densities and 
curvatures to be small enough to justify working within a semi-classical analysis; (ii) we 
compute brane-bulk couplings to leading order in a derivative expansion (making the domi- 
nant players the brane tensions and brane-localized fluxes); (m) we seek solutions that are 
axially symmetric in the two dimensions transverse to the branes and whose on-brane geome- 
tries are maximally symmetric; and (iv) we linearize the brane properties about the choices 
that source simple rugby-ball geometries (which have flat on-brane geometries despite having 
nonzero tensions). 

The last two of these assumptions deserve some motivation. The linearization about 
rugby ball geometries || [l2|, [l4|, [l^] is made in order to allow the search of their immediate 
neighborhoods in field space to be systematic; the linearity of the equations allows the con- 
struction of their most general solutions. We do not believe that the qualitative features of 
our results (like the existence of very large volume solutions, and the suppression of on-brane 
curvatures) depend strongly on this assumption. 

By contrast, at first blush the assumption of maximal symmetry might seem more re- 
strictive, since maximal symmetry is known not to be possible for a majority of brane config- 
urations and the general situation is expected to be time-dependent Jl7|]. We employ a trick 
to explore such configurations: we stabilize the time-dependent runaway by turning on an 
external current that couples to the system's low-energy moduli. In this way we can explore 
the potential energy cost that drives these runaway solutions, at least at the low energies of 
main interest. 

The supergravity of interest has a one-parameter flat direction, labeled by a particular 
combination of the 6D dilaton and the breathing mode of the extra-dimensional metric. Our 
main calculation interest for this theory is in the potential energy generated for this flat 
direction by the back-reaction of the bulk-brane couplings; and in the related change of shape 
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of the extra-dimensional geometry. We find that these display the following noteworthy 
features: 

• Volume stabilization: Any non-derivative coupling of the branes to the bulk dilaton, 
<f>, breaks the classical scaling symmetry of the bulk field equations, and so lifts the 
degeneracy of the classical zero mode. Because the bulk volume depends exponentially 
on the canonically normalized dilaton, </>, we find that a mild hierarchy in the brane- 
bulk coupling parameters can easily generate an exponentially large extra-dimensional 
volume. The exception to this is if the branes also couple only to exponentials of </>, as 
is in particular often true for D-branes. 

• Suppressed on-brane curvature: A remarkable feature of rugby ball geometries is that 
their on-brane directions are flat despite the presence of large brane tensions. We find 
that perturbations about these geometries can - but need not - share this feature, 
having on-brane curvatures that are parametrically small compared with the generic 
size of the on-brane tensions. In particular, a mechanism for achieving such solutions 
arises for some types of dilaton-brane couplings since the dilaton can be driven to roll 
out to large fields along the flat direction to find places where the on-brane tension and 
curvature are the smallest. 

• Relevance of on-brane fluxes: Flux quantization within the bulk provides a strong con- 
straint on flux-stabilized rugby-ball geometries, and in the simplest examples gives 
rugby-ball perturbations whose on-brane curvatures are not suppressed. An important 
part of our ability to find other solutions with lower curvature is our inclusion of brane- 
localized flux, corresponding to a magnetic coupling of the branes to the geometry- 
stabilizing fluxes. In this way our calculations bear out the earlier expectations of 
ref. §. 

• Quantum corrections: In section [D] we provide a preliminary estimate of the size of 
quantum corrections for a particularly promising toy model, with both exponentially 
large volumes and a suppressed on-brane curvature. Our estimates indicate these prop- 
erties need not be destabilized by quantum effects on the brane or in the bulk. They do 
not do so because it is the value taken by the dilaton along the classical flat direction 
itself that plays the role of the loop-counting parameter (similar to what happens for 
string vacua), and this constrains how quantum effects can alter the dynamics that 
determines what this value is. In particular, we find that provided the brane couplings 
are arranged to lie within the regime of weak coupling, the conclusions of the classical 
analysis are protected from loop effects. 

• Bulk distortion: Although the rugby ball solutions themselves are simple we find that 
the nearby geometries are more generic, including warping and nontrivial dilaton profiles 
across the extra dimensions. 
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• Modulus-matter couplings: Having the dilaton couple to branes only through an overall 
prefactor implies a universal coupling to ordinary matter, if this resides on a brane. 
This is likely to have interesting phenomenological implications if the moduli can be 
arranged to be light enough to mediate macroscopic forces. It is noteworthy that these 
couplings can have the form required to profit from a 'chameleon' mechanism [p^j. 

We regard these properties to be new examples of how low-energy brane dynamics can 
change how one thinks about technical naturalness and hierarchies of scale. In particular, the 
natural generation of exponentially large volumes in the 6D model explored here fills in a key 
missing step in efforts to use large volumes to solve the gauge hierarchy problem. 

Although none of the solutions explored here have on-brane curvatures that are low 
enough to describe the Dark Energy density, the existence in some models of a mechanism for 
robust parametric suppression of the on-brane curvature is very suggestive. We regard this as 
encouragement to continue to explore this direction for new approaches to the cosmological 
constant problem. 
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A. Flux quantization with brane fluxes 

To see how to interpret the parameter rewrite the brane flux term as a regularized 6D 
integral weighted by a scalar function s(p) whose support is nonzero only in a short interval 
\p — Pb\ < £ away from the brane, and is normalized so that J d 2 x yfg~2 s = 1. That is, 





Then the 5A.Q Maxwell equation becomes 





which integrates to give 




(A.3) 
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This is the bulk solution found in the text away from the brane, where s = 0. 

Imagine now integrating this to obtain Ag(p) in the vicinity of the brane at pb = 0, using 
for s a simple step function: s = l/(ire 2 ) for p < e and s = for p > e. Assuming W ~ Wb is 
approximately constant and e B ~ p for p < e, we set A$(0) = and integrate tof find Ag(e), 



Me) = ~2 



1 2 

2 P 



+ Q [ dppe* 
o Jo 



(A.4) 



and so as long as diverges less fast than p 2 at the brane we have 

UmA e (e) = p-. 

e->0 ITT 



(A.5) 



The junction condition for A' g at p = e can also be seen by subtracting the solution, 
eq. ( |A.3| ) evaluated at p < e — where s = l/(7re 2 ) — from the same solution evaluated at 
p > e — where s = 0. Since the RHS is the same in both cases we get the following jump 
discontinuity across p = e: 

p=e+ $ 5 



t2 ' 



(A.6) 

p=£- TTE* 

This can be related to the derivative of the brane action with respect to Aq by rewriting 



eq. ([A.lD as 



$6 J d 6 x s T pd = ^ J d 4 x A e (e) , (A.7) 



and so (keeping in mind the relative sign between the tension and flux terms) 



e~ B A' ft 



p=£+ 
p=£- 



1 



2vr \dAe 



(A. 



as stated in ref. 

B. Alternative currents 

In this section, we check that the details of the current are not important for the stabilization 
of (p. Define for comparison purposes the current 



Sj 



d 6 x 



(B.l) 



This choice keeps the scale invariance of the bulk action intact. Compared to the current used 
in the main body, the changes to the linearized equations of motion arise only in the <j) and 
SB equations. Here we write only the contributions to these equations due to the current: 



(sinxScj)') 
(sin 2 x 



SB')' 



sin 2 x 



r^JL^e^ sin a; 
k 2 Je w , 



(B.2) 
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where the SB term was present before but lacked the factor e vo . By contrast the 5(f) contri- 
bution given above didn't exist for the current used in the main text. 

The resulting change to the perturbations that is a consequence of the current only is 



(5(f>) j = k 1 JL 2 e^° In Isinx 
(SB) j = K 2 JL 2 e^° 



-%2(x) + - (1 + X COt x) 



(B.3) 



The resulting changes in the matching conditions to the brane are as follows: 



Q 4ira 



-K 2 JL 2 e^ 
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5T N + 5T S - Q$ N - Q$ s + - (T' N + T' s ) 



(B.4) 



Relative to the main text, the current appearing in the first of these equations is new, and in 
the second equation it is half the size as found in the main text (apart from the trivial scaling 
by e^ throughout). The resulting 4D curvature is unchanged because it is the combination 



1}r = 2 - 2 K 2 JL 2 e^ = ^-(t' n + T' s ) - An 2 JL 2 e^ 



(B.5) 



Here the current contribution is twice the result of the main text, so with the current being 
only half as large for a given ip , the curvature remains unchanged. 

Using the corresponding current in the 4 dimensional theory — i.e. using \J — g j — yields 
in general a different Einstein- frame effective potential. However, Vef^*) and Vj^- (</?*) agree. 
Since the coefficient of the kinetic term is unchanged, neither is the mass of the dilaton and 
the cosmological constant. This shows that we can extract the properties at the stationary 
point reliably, even though the shape of the potential away from this point can depend on the 
detailed definition of the current that is used. This reflects a general property: the detailed 
form of a scalar potential can be varied (as always) by performing a field redefinition, though 
any dependence on the field variables used ultimately drops from any physical prediction. 



C. Linearization around the rugby ball 

This appendix computes the linearization of the field equations about the rugby ball solutions, 
with 



B 



e Bo (l + SB) = e~^l 2 aL sin ( £ ) (1 + SB) 



W = 5W and <f> = (po + S<f> (C.l) 
T pd = (Q + SQ)e 4l+B - m = Qe^ +Bo (l + ^ + SB + 5(f> - ASW^j . 



- 34 - 



Using these in the d'Alembertian for the dilaton gives 
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(C.2) 



where we use <9p</> = d p (5cp) to allow the use of the background metric. Similarly, 
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where the second line uses the rugby ball condition for the background value of Q. With 
these the dilaton equation becomes 
dp[sm(p/L)dp{6<t>)] 



sin (p/L) 

which is the form used in the main text. 

Similarly, the linearization of the Einstein equation, eq. (p 
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Since J is perturbatively small, the background metric can be used to simplify the current 
term in the action, 



<*(\/=S J) = 
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ensuring that J appears as a new contribution k 2 J to the Einstein equations, eq. (2.7). 
Finally, the derivative terms for 5B become 
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d p [ S m 2 (p/L)d p (5B)} 
sin 2 (p/L) 



(C.7) 



Putting this all together yields 
d p [sin 2 (p/L) d p (5B)] _J_ 

sin 2 (p/L) L 2 
dp [sin 2 (p/L) d p (SB)] _ l_ 

sin 2 (p/L) L 2 



■ + 



• - ( • ; (>t]V +K 2 JL 2 e~ lp0 



5Q 
Q 

5Q 
Q 



and 



R 



2W 1 

U + L COt 



d p W + djW 



65W + K 2 JL 2 e~' fi0 
2e^° fSQ 



+ 



L 2 \Q 



4 

— cot 

Li 



d 2 ,W. 
p 



2k 2 J, 



d p w 



(C.8) 
(C.9) 



which are the equations solved in the main body. 
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D. Some useful integrals 



This appendix evaluates the integrals Mi and Hi encountered in the main text. 
Evaluation of M 

This section evaluates the constant encountered in the ip\ perturbation. The integrals of 
interest are 



Mi 
M 2 
and M 



(x) = r 

Jo 

(x) = f 
J 



dy sin y In 



cos y 



smy 



dy — 



MM 



sin y 



PIT 

/ dx sinx M 2 (x) . 
Jo 



(D.l) 



First of all, notice that the logarithm in the first line is antisymmetric under y — > tt — y, 
while sin 2 y is symmetric. This means that the first integral integrates to if x = it: that is, 
Mi{ir) = 0. These observations justify the following manipulations: 



Mi(-k - x) = Mi{ix) 



dz sin zln 



= — dy sin 2 (7r — y) In 



1 — cos z 



sin z 
1 — cos(-7r — y) 



sin(7r — y) 



J o 



dy sin y In 



1 — cos y 



smy 



The same manipulations applied to A^2 then give: 

M 2 (ir - x) = M 2 (tt) - 

J n~x 



(D.2) 



dz 



Mi(z) 
sin 2 z 



= M 2 {tt) - / dy 
Jo 

= M 2 (tt) - dy 
Jo 



sm 2 (vr - y) 
Mi(y) 



sin 2 y 



= A^ 2 (vr) - M 2 {x) . 



(D.3) 



Numerical evaluation of M 2 {tt) is complicated by the weak convergence of the integral 
near tt. It can be evaluated more efficiently by using the above expressions to relate it to 
M 2 {it/2). That is, numerical integration gives M 2 (tt/2) = —0.5 to within the numerical 
(Maple 11) precision. Using this, we find 



M 2 (ir) = 2M 2 (n/2) = -1. 



(D.4) 
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Hence M.2 satisfies 

M 2 (tt - x) = -1 - M 2 (x) . (D.5) 

To evaluate M, use 

fir/2 



rn/Z rn 

M. = / dx sinxA^Oc) + / dx sinxA^OzO 

Jo Jtt/2 
rn/2 r-n/2 

= dx sinxA^x) + / dx sin(7r — x) (vr — x) 

Jo Jo 

rn/2 r-n/2 

= dx sinxA^x) + / dx sinx — 1 — M. 2 {x) 

Jo Jo L J 

r-Tr/2 

= — dx sinx = — 1 . (D.6) 
Jo 



Evaluation of % 

Recall the definitions, 



rx 

%\{x) = dy sin 2 y In | sin y \ 
Jo 

n, ( , f x A ^i(y) 

U 2 {x) = / dy . 2 
Jo sin y 

W = Tdy s±nyU 2 {y). (D.7) 
Jo 

In this case numerical evaluation gives (Maple 11): 

Hi(7r) = ^(l-ln4). (D.8) 

Similar numerical integration to evaluate 'H 2 ((f>) is complicated by the apparent singularity 
at the endpoints caused by the factors of 1/ sin 2 y in the integrand. These can be dealt with 
by repeating the arguments of the previous section, which in this case give 



7T)- / . 

J TT — X 



Hi(n — x) = y,i(ir) — I dy sin y In |siny 

t — : 
rx 

= "Hi(vr) — dy sin 2 y ln|siny| 
Jo 

= Hi(ir) -Ui{x). (D.9) 
Next consider the following symmetry properties of H 2 : 

rn/2+x 



fit \ /vr\ r< z+x 



n/2 sin 2 z 



n/2 sin 2 (7r-y) 

d, (d,o, 

v 2 / sm z 2/ 
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and simplify using 



* /2 , um r'\ %M r' 2 ~ x %m 
d y ~~t:. = / & y ~ / d y —2— 



tt/2-x sm"y Jo sury J sm~ y 

= „ 2 g)_„ 2 (|_ x ), (D . n) 



to get 



«.(5 + »)=%(5-»)+«.W r 

V 2 / V 2 J Jn/2-x sm z y 

= « 2 (|-x) +Wi(tt) cot(|-x) . (D.12) 

The evaluation of "H now proceeds, with 

/>7T 

H = dx sinxT^x) 
Jo 

= J dx jsinxT^Oz) + sin ^ + x^j H2 (J^ + x^j j 

= y dx |sinx^ 2 (^) + sin ^ - x^j H2 - x^j + Hi(-k) cot ^ - x) j | 

= 2 / dx sinx^2(^) + "Hi(tt) / dx cos ^— — xj 

/•tt/2 

= 2/ dx sinx^ 2 (a;) +Hi(tt) . (D.13) 

JO 



This can now be integrated numerically without problems, giving (to ten decimal places) a 
result consistent with % = — 2 + In 4. 
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